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Computer tomography method using redundant measured values 



The invention relates to a computer tomography method in which the 
measured values used to reconstruct a CT image are only ones that lie within a reconstruction 
window, and it relates in particular to exact computer tomography methods. The invention 
also relates to a computer tomograph for performing the said method and to a computer 
5 program for controlling the computer tomograph. 

In known methods of the above-mentioned kind, a radiation source moves 
along a helical or circular trajectory relative to an examination region in which the object to 
be examined is situated. As it does so, a detector unit acquires measured values from which a 
CT image, such as, for example, an image of the distribution of absorption in the examination 

10 region, can be reconstructed. Many reconstruction methods, and particularly the exact 
reconstruction methods, are able to use only a proportion of the measured values. The 
proportion of the measured values that can be used for reconstruction may for example be 
defined by a range of solid angles through which the radiation source moves while a point in 
the examination region that is to be reconstructed is being irradiated, or by a given region on 

1 5 the detector unit. The region in which measured values that are used for reconstruction are 
situated is called a reconstruction window. Measured values outside the reconstruction 
window are ignored by methods of the kind mentioned and are termed redundant measured 
values. A reconstruction method of this kind has been described by, for example, Katsevich 
in "Analysis of an Exact Inversion Algorithm for Spiral Cone-Beam CT", Physics Medicine 

20 and Biology, vol. 47, pp. 2583-2597 (El). In this case, the radiation source moves along a 
spiral, i.e. helical, trajectory relative to the examination region and the measured values that 
are used are only the ones that lie within what is termed the PI window, which will be 
explained in detail below. Things that are disadvantageous about such methods are that they 
have a poor signal-to-noise ratio and that, if the object moves, they produce motion artifacts, 

25 which are a nuisance. 

It is therefore an object of the present invention to specify a computer 
tomography method of the kind stated in the opening paragraph in which the sighal-to-noise 
ratio is improved and the motion artifacts are reduced. 
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This object is achieved in accordance with the invention by a computer 
tomography method having the following steps: 

a) generation, by a radiation source, of a conical beam of rays that passes through 

an examination region and an object situated therein, 
5 b) production of a relative movement between the radiation source on the one 

hand and the examination region on the other hand, which movement comprises at least a 
rotary movement about an axis of rotation and is in particular in the form of a helix or circle, 

c) acquisition, by a detector unit and during the relative movement, of measured 
values that depend on the intensity in the beam of rays on the farther side of the examination 

10 region, 

d) determination, with the help of redundant measured values, of a 
complementary measured value for each of at least some of the measured values that were 
acquired in step c) and that lie within a reconstruction window, the rays associated with the 
given measured value and the complementary measured value belonging to it being oriented 

15 in opposite directions to one another, 

e) replacement of each measured value for which a complementary measured 
value was determined in step d) by a sum comprising the measured value, having been 
weighted, and the complementary measured value, having been weighted, 

f) reconstruction of a CT image of the examination region from the measured 
20 values lying within the reconstruction window. 

In contrast to known computer tomography methods of the kind defined in the 
opening paragraph, redundant measured values are used to determine complementary 
measured values that, when added to their respective associated measured values, replace the 
latter, which means that redundant measured values too contribute to the reconstruction that 

25 follows. The invention is based on the finding that, in computer tomography methods in 

which the reconstruction is not confined to only a proportion of the measured values, motion 
artifacts are reduced and the signal-to-noise ratio is improved. Taking account of the 
redundant measured values therefore results not only in an improvement in the signal-to- 
noise ratio but also in a reduction in the motion artifacts. 

30 In claim 2, if the complementary measured value has not been acquired in step 

c) then it is determined by means of John's equation. This produces CT images of particularly 
good quality. 

John's equation is a known partial differential equation that, among other 
things, allows further measured values to be determined for measured values that have been 
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acquired, which further measured values are consistent with those that have been acquired. 

The term "consistent" can be explained as follows. The object in the 
examination region is not clearly defined by the set of measured values acquired in step c) 
because the measured values are affected by noise or also, for example, because too few 
5 measured values have been acquired. For the set of measured values, there are therefore a 
number of possible objects in the examination region that fit the measured values. A 
complementary measured value is consistent with the measured values that have been 
acquired when it can be produced by the passage of a ray through one of the possible objects. 
A consistency condition with which a complementary measured value of this kind can be 

10 determined is for example John's equation, which is known. 

In claim 3, when the weighted complementary measured values and the 
respective weighted measured values belonging to them are added, a complementary 
measured value that has been obtained is given a heavier weighting if it was acquired in step 
c) rather than having been determined by, for example, John's equation. The greater account 

15 that is thereby taken of values that have actually been measured gives images of an improved 
quality. 

Claims 4 and 5 describe a preferred exact computer tomography method that 
produces CT images of good quality for helical trajectories. 

A computer tomograph for performing the method according to the invention 
20 is described in claim 6. Claim 7 defines a computer program for controlling the computer 
tomograph claimed in claim 6. 

These and other aspects of the invention are apparent from and will be 
elucidated with reference to the embodiments described hereinafter. 



25 

In the drawings: 

Fig. 1 is a schematic view of a special embodiment of computer tomograph 
with which a special embodiment of the method according to the invention can be performed. 

Fig. 2 is a flow chart of the method according to the invention. 
30 Fig.3 shows a PI straight line and a PI interval for a point in the examination 

region. 

Fig. 4 shows a PI straight line and the PI interval for a point in the examination 
region, when projected onto a plane perpendicular to the axis of rotation. 

Fig. 5 is a schematic view of a planar detector having a PI window. 
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Fig. 6 is a schematic view in cross-section of a helical trajectory and of an 
examination region, showing a direct ray and a complementary ray. 

Fig. 7 is a flow chart for an exact reconstruction of a CT image. 

Fig. 8 is a schematic perspective view of parallel rays in different positions. 
5 Fig. 9 is a schematic perspective view of the helical trajectory showing a k 

plane and a k line. 

Fig. 10 is a schematic perspective view of the helical trajectory showing, in 
parallel planes, fans of rays formed by rebinning. 

10 

The computer tomograph shown in Fig. 1 comprises a gantry 1 which is able 
to rotate about an axis of rotation 14 extending parallel to the z direction of the coordinate 
system shown in Fig. 1 . For this purpose, the gantry 1 is driven by a motor 2 at an angular 
velocity that is preferably constant, but adjustable. Fastened to the gantry 1 is a radiation 
15 source S such as an X-ray generator, for example. This is provided with a collimator 

arrangement 3 that extracts from the radiation produced by the radiation source S a conical 
beam of rays 4, i.e. a beam of rays that is of a finite extent other than zero both in the z 
direction and also in a direction perpendicular thereto (i.e. in a plane perpendicular to the axis 
of rotation). 

20 The beam of rays 4 passes through a cylindrical examination region 13 in 

which an object, e.g. a patient on a patient presentation table (neither of which is shown) or 
even a technical object, may be situated. Having passed through the examination region 13, 
the beam of rays 4 impinges on a detector unit 16 fastened to the gantry 1, which detector 
unit 16 has a face that comprises a plurality of detector elements that, in this embodiment, are 

25 arranged in the form of a matrix in rows and columns. The columns of detector elements 

extend parallel to the axis of rotation 14. The rows of detector elements are situated in planes 
perpendicular to the axis of rotation and, in this embodiment, are situated on an arc about the 
radiation source S (thus forming a focus-centered detector unit). In other embodiments 
however, they may be differently configured, e.g. to describe an arc about the axis of rotation 

30 14 or in a straight line. In any position of the radiation source, each detector element on 

which the beam of rays 4 impinges gives a measured value for one ray from the beam of rays 
4. 

The included angle of the beam of rays 4, which is designated ow, determines 
the diameter of the object cylinder within which the object being examined is situated when 
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the measured values are being acquired. The included angle is defined in this case as the 
angle that is made with a plane defined by the radiation source S and the axis of rotation 14 
by a ray that, in a plane perpendicular to the axis of rotation 14, is situated at the edge of the 
beam of rays 4. The examination region 13, i.e. the object or the patient presentation table, 
5 can be displaced by means of a motor 5 parallel to the axis of rotation 14, i.e. to the z axis. 
However, to give an equivalent action, the gantry could equally well be displaced in this 
direction. When the object is a technical one rather than a patient, it can be rotated in the 
course of an examination while the radiation source S and the detector unit 16 remain 
stationary. 

10 If the motors 2 and 5 run simultaneously, the radiation source S and the 

detector unit 16 describe a helical trajectory relative to the examination region 13. If on the 
other hand the motor 5 for advance in the direction of the axis of rotation 14 does not run and 
the motor 2 causes the gentry to rotate, a circular trajectory is produced for the radiation 
source S and the detector unit 16 relative to the examination region 13. In the present 

15 embodiment, only the helical trajectory will be considered. The method according to the 
invention could, however, also be applied where the trajectory was circular. 

The measured values acquired by the detector unit 16 are fed to a 
reconstructing computer 10 that is connected to the detector unit 16 by, for example, a data 
transmission system (not shown) that operates without electrical contacts. The reconstructing 

20 computer 10 reconstructs the CT image and reproduces it, on a monitor 1 1, for example. The 
two motors 2 and 5, the reconstructing computer 10, the radiation source S and the transfer of 
the measured values from the detector unit 16 to the reconstructing computer 10 are 
controlled by a control unit 7. 

Fig. 2 shows the flow of a method of measurement and reconstruction that can 

25 be performed with the computer tomograph shown in Fig. 1 . 

After initialization in step 101, the gantry rotates at an angular velocity that is 
constant in the present embodiment. The velocity may however equally well vary, e.g. as a 
function of time or of the position of the radiation source. 

In step 103, the examination region, i.e. the object or the patient presentation 

30 table, is displaced parallel to the axis of rotation and the radiation from the radiation source S 
is switched on, thus enabling the detector unit 16 to detect the radiation in a plurality of 
angular positions. In the course of this, the radiation source S moves relative to the 
examination region along a helical trajectory of which the parameters can be defined by 
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OJcos^ 
y(A.)= RsmX 



(1) 



Here, 2ith is the feed of the table per rotation, X is the angular position of the 
radiation source relative to a reference angular position that is any desired angular position 
but fixed, and R is the radius of the helical trajectory 17. 

In the present embodiment, the advance per rotation is selected to be such that, 
from any point x in the examination region, the radiation source S is visible over an angular 
range of at least 180°. The picking up of measured values with a table advance of this kind is 
referred to as PI acquisition 

PI acquisition will be explained in detail in what follows. For this purpose, 
there are shown in Fig. 3 the helical trajectory 17 along which the radiation source moves 
relative to the point x in the examination region, and the section /pi(x) of the helix 17 that is 
cut off by a PI straight line 37. Fig. 4 shows the helix 17 from Fig. 3 projected onto a plane 
oriented perpendicularly to the axis of rotation. The PI straight line 37 is that line which cuts 
the helix at two points and the point x, with the section ipi(x) of the helix that is cut off by the 
straight line extending through an angle of less than 2n. In PI acquisition, rays that start from 
the radiation source S while the latter is on section /pi(x) of the helix pass through the point x. 

In the context of the present embodiment, what is used below is a method of 
reconstruction that employs only measured values whose associated rays originate from 
section ipi(x) of the helix, i.e. the measured values used are only ones that lie in what is called 
the PI interval. It is known that, by considering the geometry, it can be shown that measured 
values whose rays start from a given section of the helix lie within a reconstruction window 
on the detector unit. In the present PI acquisition, the reconstruction window is the PI 
window. The PI window 25 is bounded on a planar detector 60 by the PI lines 2 1 , 23 that are 
defined by the following equations: 



Here, u and v are coordinates on the planar detector 60 in the coordinate 
system shown in Fig. 5 . For reasons of clarity, the coordinate system has been shown below 




(1) 



and 




(2) 
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the planar detector 60. However, the origin of the coordinate system is situated at the center 
of the detector. The unit vectors of the coordinate system are l u = (- sinX, cosX, 0) and l v = (0, 
0, 1). 

The planar detector 60 is an imaginary detector that contains the axis of 
5 rotation 14 and is oriented perpendicularly to that ray that, starting from the position of the 
radiation source at the time, strikes the axis of rotation 14 perpendicularly. Measured values 
that are detected by the actual detector unit 16 can be projected onto the planar detector 60 
along their corresponding rays. 

Computer tomography methods of the kind defined in the opening paragraph 

10 use only measured values that lie within a reconstruction window. The rest of the measured 
values are not used and this leads to a poor signal-to-noise ratio and motion artifacts. In the 
following steps, a complementary measured value gc(ko,UQ,vo) is determined, for each 
measured value g(Xo,w<j>vo) that lies within the reconstruction window, with the help of 
redundant data, i.e. measured values that lie outside the reconstruction window, as a result of 

15 which the signal-to-noise ratio is improved and the motion artifacts reduced. The ray 
associated with the measured value g(ko 9 uo 9 vo) is referred to below as the direct ray. 

In other embodiments, a complementary measured value may be determined 
for each of only a proportion of the acquired measured values that lie within the 
reconstruction window. In particular, it would be possible for a complementary measured 

20 value to be determined only for those measured values that are situated in the edge regions of 
the reconstruction window. 

A complementary measured value is defined by the path followed by the 
complementary ray belonging to it. A complementary ray follows a path parallel to the direct 
ray belonging to it but in the opposite direction thereto. A consistency condition can be used 

25 to determine the complementary measured value associated with the complementary ray. The 
consistency condition that is used in the present embodiment to determine the complementary 
measured value associated with the complementary ray is the equation generally referred to 
as John ! s equation that was published in "Partial Differential Equations", F. John, Applied 
Mathematical Sciences, Springer- Verlag, 1971 and that has been adapted for helical 

30 trajectories in US 6,292,526 B 1 . 

In accordance with the invention, the method used to determine a 
complementary measured value may be any method, and in particular any consistency 
condition, that, by taking account of redundant measured values, allows a complementary 



WO 2004/111945 



8 



PCT/IB2004/050845 



10 



15 



20 



25 



(3) 



measured value to be determined, at least approximately, for a complementary ray that 
follows a path parallel to the direct ray but in the opposite direction thereto. 

To determine the particular complementary measured value g c (^o»w<>,vo), ^ 
step 105 an approximating measured value g(Xi,Wi,vi), whose associated approximating ray 
3 1 follows approximately the same path as the complementary ray (see Fig. 6), is first 
determined from the acquired measured values. This approximating ray starts from the 
angular position XI and impinges on the planar detector 60 at a point (u\, vj). 

The angular position of the approximating ray is defined by 

fA, 0 +7i-2p where v 0 >0 
[A, 0 -7t-2p where v 0 <0 

Here, P is the fan angle of the direct ray, i.e. the angle that the direct ray makes 
with a ray that intersects the planar detector 60 in Fig. 5 in the center. The angular position Xi 
is the angular position at which a direct ray that starts from the position y(Xo) of the radiation 
source and impinges on the planar detector at a point (u 0> v 0 ), first touches the helix 17 if the 
ray is displaced parallel to the axis of rotation 14. 

The equation that governs the u coordinate of the approximating ray on the 
planar detector 60 is: 

u\ = -w 0 (4) 
The v coordinate is determined as follows. By X = X } and u = u\ is defined a 
fan of rays that starts from the angular position X\ on the helix 17 and impinges on all those 
detector elements for which u - u\ is true. The respective point of intersection with the direct 
ray 33 is determined for each ray in this fan of rays. The approximating ray 3 1 is then that ray 
in the fan of rays whose point of intersection is closest to the axis of rotation 14. This 
approximating ray impinges on the planar detector at a point (i*i,vi), where 



v, = 





f 

A(2P-7l) 

V 


R 


V«> 2 +* 2 


f 

/j(2p+7t) 


R 





R 



2cosP 

1 

2cosP R 



-2v 0 cosP 



-2v 0 cosp 



where v 0 >0 



where v 0 <0 



(5) 



In step 107, the path followed by the complementary ray is determined. As has 
already been explained above, the complementary ray follows a path parallel, but in the 
opposite direction, to the direct ray 33. The direct ray 33 impinges on the surface 35 of the 
cylinder that is defined by the helical trajectory 17 at a point 



WO 2004/111945 



9 



PCT7IB2004/050845 



z=z 0 =hX 0 +v 0 . (6) 

<sju 2 +R 2 

During the acquisition of the approximating measured value, the radiation 
source S was situated at a point 

h(X 0 + n-2$) where v 0 >0 

lh(X Q -n~2$) where v 0 <0 j 

To allow a complementary ray to be obtained, the position z\ of the radiation 
source has to be "shifted" to the position z 0 . It must therefore be shifted by 



Aq=z 1 -z 0 = 



A(+*-2p)-v 0 ^ CQS P W here v 0 >0 
/,(_^.2P)-v 0 ^ C ° S|3 where v 0 <0 



The position of the radiation source for the complementary ray is thus that 
position of the radiation source which is defined by the angular position X\ and has been 
10 shifted along the axis of rotation by A<;. 

The complementary ray that starts from that radiation source position of the 
approximating ray which has been shifted by Ac; must impinge on the surface of the cylinder 
at the z position from which the associated direct ray starts. This produces the following 
defining equation for the coordinate vi of the complementary ray: 

15 hX 0 =hX l+Vi **pl (9) 

By inserting equating (3), this gives 



h(2Q-n)-*-^ where v 0 >0 

v H ' 2i?cosp 0 

lu 2 + R 2 

A(2p+7i)^ — where v 0 <0 

2R cosp 



(10) 



The complementary ray is therefore found from the approximating ray by 
shifting the radiation source along the axis of rotation by Aq and by shifting the point of 
20 impingement of the ray on the detector by 



R 2 



Av=v 1 -v, = -Aq — ^2*- 01) 



The u coordinate of the complementary ray is equal to the corresponding 
coordinate of the approximating ray. 
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Once the path followed by the complementary ray complementary to a direct 
ray is known, the complementary measured value associated with the complementary ray can 
be determined in step 109. If the complementary measured value was acquired, the measured 
value acquired can be used as a complementary measured value. This is not generally the 
5 case however, because the radiation source position of the complementary ray is usually not 
situated on the helix 17. 

From "Improved 2D rebinning of helical cone-beam CT data using John's 
equation", Defrise, Noo, Kudo, M10-74 in 2002 Nuclear Science Symposium Conference 
Record, guest editor Scott Metzler, 10-16 November, Norfolk, Virginia, USA, ISBN 0-7803- 
10 7637-4, it is known that, with the help of John's equation, a complementary measured value 
that is consistent to a good approximation with the measured values can be calculated by 
adding a correcting measured value to the approximating measured value. The correcting 
measured value Ag(X u uuv\) is given by the following equation: 




15 In accordance with the invention, all or only a proportion of the measured 

values acquired can be used to determine the correcting measured value from equation (12). 
In the present embodiment, the measured values used are only ones whose associated rays 
start from radiation source positions y(X) that are located on the helix 17 in the vicinity of the 
radiation source position y(ki) of the approximating ray. What that means is that the 

20 measured values used are only ones whose associated radiation source positions are situated 
in the range y(X) where X\ - Ax < X < X\ + Ax- In the present case, Ax is selected to be such that 
ten different radiation source positions y(X) are covered. In other embodiments a different 
number of radiation source positions, and hence of measured values, can be used. 

The range Ax is selected to be at least sufficiently wide for redundant measured 

25 values to contribute to at least some complementary measured values. However, in the 
majority of cases this condition will already have been met due to the fact that the 
approximating measured value is itself redundant. 

To determine the correcting measured value, all the measured values whose 
associated rays start from radiation source positions that are within the above range are 

30 derived from equation (12) by their angular positions X and their v coordinate, to allow a first 
dataset of partial derivatives to be derived from equation (12). These derivations, and the 
following ones too, may be performed by, for example, the finite difference method. Also, 
the original measured values, i.e. the non-derived ones, whose associated rays start from 
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radiation source positions that are within the above range are twice derived as partial 
derivatives, for v, to allow a second dataset of derivatives to be obtained. 

The measured values in the second dataset of derivatives are then multiplied 
by a weighting factor that follows a square law in u and is linear in v. In the present 
5 embodiment, this weighting factor is equal to 




The next step is to add the measured values in the first and second datasets of 
derivatives to give a set of measured values, and those resulting values whose associated rays 
start from the angular position X\ and impinge on row vi of the detector unit and on columns 

10 whose coordinates u are equal to or smaller than u\ 9 are added up. 

Finally, the resulting sum is multiplied by a factor that is proportional to A<; 
and inversely proportional to R. The sum is preferably multiplied by Ac/R. The sum 
multiplied by the correcting factor is the correcting measured value for the approximating 
measured value g(Xi»tfi»vi), or for the measured value g(Ao,wo,vo) of the direct ray. 

15 The correcting measured value having been calculated, the complementary 

measured value g^o^vo) to the measured value g(^o 9 u 09 vo) can be determined from 

gc(Xo,wo,v 0 ) = gfruuuvi) + AgiX U UuVy) (14) 
by adding the correcting measured value to the approximating measured value. 

Once steps 105 to 109 have been performed for all the measured values 

20 gQw>uo 9 vo) that lie within the PI window, two sets of measured values exist, i.e. the dataset 
acquired in step 103 and the dataset complementary thereto. 

In step 111, each measured value for which a complementary measured value 
has been determined is replaced by a sum comprising the complementary measured value 
after it has been weighted and the associated measured value after it has been weighted. For 

25 this purpose, the given measured value is added to the associated complementary measured 
value, each value having been multiplied by a weighting factor prior to the addition. If the 
complementary measured value was acquired, then it and the associated measured value are 
weighted by equal amounts, and in particular are each multiplied by 0.5. If the 
complementary measured value was determined from a consistency condition, such as John's 

30 equation for example, then it is multiplied by a smaller weighting factor than the associated 
measured value. In this way, the measured value may be multiplied by 0.9 or 0.8 for 
example, whereas the complementary measured value is multiplied by 0.1 or 0.2. In the 
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present embodiment, the sum of a measured value/complementary measured value pairing is 
equal to 1. 

In the next step 1 13, a CT image, which in this embodiment is an image of the 
distribution of absorption in the examination region, is reconstructed by an exact method 
5 from the measured values lying within the PI window that were determined in the previous 
steps, i.e. in particular from the measured values that have been obtained by adding weighted 
complementary measured values and associated measured values. The individual steps in the 
reconstruction are shown in Fig. 7. 

In other embodiments, a complementary measured value could also have been 
10 determined for only a proportion of the measured values that lie within a reconstruction 

window. If this were the case, then what would be used for the reconstruction would be, on 
the one hand, the measured values lying within the reconstruction window that had been 
obtained by adding a weighted complementary measured value and a weighted measured 
value or, on the other hand, if no complementary measured values had been determined for 
15 some of the measured values lying within the reconstruction window, the measured value that 
was acquired in the particular case. 

To allow the exact reconstruction to be understood, the following equation 
from El will first be cited: 



(15) 

9=X 



± Jdr 1 8 

/pi(x) 

20 This equation describes an exact reconstruction of absorption by back- 

projection of the measured values. In it,f(x) is the spatial distribution of absorption in the 
examination region at point x. 

The measured value DJ(y,0) can be defined by the line integral: 

D / (y(X),0)=Jd//(y+/e) (16) 
o 

25 Here, the unit vector © gives the direction of the ray belonging to the 

measured value. 

In step 201, the measured values are derived as partial derivatives from 
equation (15) for q, i.e. for the angular position of the radiation source at the point q = A. It 
should be borne in mind in this case that only y depends on q and not 0, which means that, 
30 for the derivation, account has to be taken of measured values from parallel rays in each case. 
Because parallel rays have the same cone angle, then, as shown in Fig. 8, in the case of the 
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focus-centered detector 16 that is used in the present case, the parallel rays 51 impinge on the 
same row 53 in the detector. In this case the cone angle of a ray is the angle that the ray 
makes with a plane perpendicular to the axis of rotation 14. For the partial derivation, the 
measured values may first be resorted. For this purpose, measured values that belong to 
5 parallel rays, i.e. to the same row 53 in the detector but to different angular positions A*,, Xb 9 Xc 
of the radiation source, are combined in each case to form a set. The measured values in each 
set are derived, numerically for example by the finite difference method, in accordance with 
the angular position of the radiation source. 

The unit vector 0 depends on the k angle y that can be defined by means of 
10 so-called k planes 55 (see Fig. 9). The k planes 55 will be explained below. 

To allow a k plane 55 to be determined, a function 



X l (X 9 X 2 )— 



mX 2 +(n-m)X 



,X<>X 2 <X+2n 



(17) 



mXHn ~ m)X \ X>X 2 >X-2n 



n 

is introduced that depends on non-negative values n and m 9 n>m 9 that are whole numbers. In 
the present embodiment, n is selected to be equal to 2 and m to 1 . Other values could 
15 however also be selected for n 9 m. Equation (15) would still remain exact; only the positions 
of the k planes 55 would change. Also defined are the vector function 



u(A,A 2 ) = 



y(^(XA 2 ))-y(X)]xba 2 )-y(X)] , ■ 

Iy(X l (X,X 2 ))-y(X)]x|y(X a )-y(X)] 1 

y(^)*yfr) . 

r—- — -i , A 2 — ^ 



(18) 



JyWxy(M' 
and the unit vector 

P(^) (19) 

20 The vector p points from the radiation source position y(X) to position x. To 

determine the k plane, a value X z e/ PI (x) is selected such that yX 9 y(X u X 9 %z) 9 y(X 2 ) and x lie in 
a plane. This plane is referred to as the k plane 55 and the line of intersection between the k 
plane 55 and the face 74 of the detector is referred to as the k line 57. In Fig. 9, the detector 
74 is bounded by two successive turns of the helical trajectory 17 and is of the same 

25 curvature as the helix 17. This detector has been used here as an example to make things 

clearer. Corresponding lines of intersection 57 can be determined for other detectors, such as 
the focus-centered or the planar detector. In Fig. 9 is shown a fan-shaped part of a k plane. 
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The edges of the fan meet at the location of the radiation source. Defining the k plane 55 in 
this way is equivalent to solving the equation 

(x-y(>0).u(Aa 2 )=0, K e/ woo (20) 
for X 2 . u is thus the normal vector of the k plane 55. To allow the vector function 0(X,x,y) to 
5 be determined, the vector 

e(X,x)=p(A,,x)xu(*,x) (21) 
is defined. With the definitions for p and e, the vector function ®(k,x,y) can then be stated as 
follows: 

©(A,, x, y) = cos y . ft(X 9 x) + sin y • e(X 9 x) (22) 
10 Because the two vectors p and e are oriented perpendicularly to u, the k angle 

Y gives the direction of the vector ® and hence the direction of a ray within a k plane. 

k planes and k lines are described in detail in El, which is hereby incorporated 

by reference. 

In step 203, the derived measured values along k lines are multiplied in 
15 equation (15) by a weighting factor that decreases as the sine of the k angle increases and in 
particular is equal to the reciprocal of the sine of the k angle, and are added up. For this 
purpose, a k line is determined for each point x in the examination region and for each 
radiation source position in which case, as explained above, a value X 2 e/ PI (x) is selected 
such that yX 9 y(Xi,X,X 2 ), y(X 2 ) and x lie in a plane, the k plane. The k line is then determined 
20 as the line of intersection between the k plane and the face of the detector. The 

multiplications by the weighting factor and the integrations or additions may be performed by 
means of, for example, a Fourier transformation. 

The derived and integrated measured values can be represented by the 
following equation: 

25 pW*<M>* )-^^ D Ay(^®(^r)U • (23) 

In this equation, p(y(X) 9 <I>(X,x)) are the derived and integrated measured 
values and ®(X 9 x) is a unit vector that points from the radiation source position y(X) towards 
point x in the examination region. 

The integration step that is missing in equation (15), i.e. the back-projection of 
30 the measured values, can then be defined by the following equation: 



m= ~^ J ) A Fkr (yW * (M)) - (24) 
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Before the back-projection, a rebinning of the measured values can take place 
in step 207. The rebinning causes the measured values to be resorted and reinterpolated as if 
they had been measured with a different radiation source (an extended radiation source, 
located on a part of a helix, that is able to emit fans of rays each parallel to one another). 
5 This will be explained in detail by reference to Fig. 1 0. In this Figure, 1 7 

denotes the helical trajectory from which the radiation source passes its rays through the 
examination region. Reference numeral 43 denotes a fan-shaped beam of rays that starts from 
the radiation source position So and whose rays extend along paths in a plane containing the 
axis of rotation 14. The conical beam of rays that is emitted by the radiation source at the 

1 0 position So can be thought of as composed of a plurality of plane fans of rays that are situated 
in planes parallel to the axis of rotation 14 and that intersect at the radiation source position 
5b- Of these fans of rays Fig. 10 shows only one, namely fan of rays 43. 

Also shown in Fig. 10 are further fans of rays 41, 42 and 44, 45 that are 
parallel to fan of rays 43 and that lie in planes parallel to one another and to the axis of 

15 rotation 14. The associated radiation source positions S. 2 , SLi and S\ 9 S 2 are occupied by the 
radiation source S respectively before and after it has reached radiation source position So. 
All the rays in the fans of rays 41 to 45 have the same angle of projection. The fans of rays 41 
to 45 define a beam of rays of a tent-like shape. 

The measured values obtained after the rebinning are then used to reconstruct 

20 the distribution of absorption in the examination region by means of a back-projection, 
carried out with equation (23) in the present embodiment. 

For this purpose, in step 209 a voxel V(x) within a presettable region (the field 
of view (FOV)) in the examination region, and a PI interval ipi(x) for this voxel, are 
determined. Then, in step 21 1, an angular position X within the interval ipi(x) is preset. In step 

25 213, a check is made to see whether a measured value whose ray passes through the center of 
the voxel V(x) is available for the angular position X. If no such ray can be found, then it is 
found at what point a central ray would have impinged on the face of the detector. The 
associated measured value is then calculated by interpolation from the measured values from 
adjacent rays. The measured value that can be assigned to the ray passing though the voxel, 

30 or the measured value obtained by interpolation, as the case may be, is multiplied in step 214 
by a weighting factor that becomes smaller as the distance between the radiation source y(X) 
and the point x to be reconstructed in the examination region increases. In the present 
embodiment, this weighting factor is, as shown in equation (24), equal to l/|x-y(X)|. In step 
215, the weighted measured value is summed to give the voxel V(x). In step 217 a check is 
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made to see whether all the angular positions X in the interval Jpi(x) have been looked at. If 
they have not, the flow chart branches back to step 211. Otherwise, a check is made in step 
219 to see whether all the voxels V(x) in the FOV have been covered. If they have not, the 
sequence continues with step 209. If on the other hand all the voxels F(x) in the FOV have 
5 been covered, then the absorption in the whole of the FOV, and hence the CT image, has 
been determined, and the exact reconstruction of the CT image comes to an end in step 221. 

As has already been mentioned above, in other embodiments the radiation 
source may also move along a circular trajectory relative to the examination region. If, in a 
case of this nature, the measured values used to reconstruct a CT image are only ones that lie 

10 within a reconstruction window, then in this case too complementary measured values can be 
determined for at least some measured values, in particular by means of John's equation. As 
was described above for the helical trajectory, a sum that replaces the measured value in the 
acquired dataset is formed from a measured value/complementary measured value pairing. 
Finally, a CT image can be reconstructed by known methods of reconstruction from the data 

1 5 that has been obtained in this way. 

A computer tomography method in which the radiation source moves relative 
to the examination region along a circular trajectory and the measured values used are solely 
ones that lie within a reconstruction window is known from, for example, "A fast and 
efficient method for sequential cone-beam CT", Koehler, Proksa, Grass, Medical Physics, 

20 vol. 28, no. 11, pp. 2318-2327. 
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LIST OF REFERENCE NUMERALS 



dmax Included angle 

X, X\ 9 X 2y X*, h)y Xc Angular positions of the radiation source on the helical trajectory 

Aq Shift in the z direction 

5 Av Shift in the v direction 

x Point in the examination region 

7pi(x) Section of helix 

S Radiation source 

&2, S-i, S 0 , S u Si Positions of radiation source 

10 1 Gantry 

2, 5 Motors 

3 Collimator arrangement 

4 Beam of rays 
7 Control unit 

15 10 Image processing computer 

1 1 Monitor 

1 3 Examination region 

14 Axis of rotation 
16 Detector unit 

20 17 Helical trajectory 

21,23 PI lines 

25 PI window 

3 1 Approximating ray 

33 Direct ray 

25 35 Surface of cylinder 

37 PI straight line 

41, 42, 43, 44, 45 Mutually parallel fans of rays 

51 Parallel rays 

53 Row on detector 

30 55 k plane 

57 k line 

60 Planar detector 

70 Beam of rays 

74 Face of detector 



